RM3-TH/99-8 



Penguin amplitudes: charming contributions 5 

M. Ciuchini a , E. Franco 6 , G. Martinelli 6 and L. Silvestrini c 

ON 

a INFN - Sezione di Roma III and Dipartimento di Fisica, 
Ph| Universita di Roma Tre, Via della Vasca Navale 84, 1-00146 Roma, Italy 

^ • 6 INFN - Sezione di Roma and Dipartimento di Fisica, 

Universita di Roma "La Sapienza", P.le A. Moro 2, 1-00185 Roma, Italy 
CN . c Technische Universitat Miinchen, Physik Department, 

D-85784 Garching, Germany 

> . 

O . Abstract 

m 

I/"") \ We briefly introduce the Wick-contraction parametrization of hadronic 

■ matrix elements and discuss some applications to B and K physics. 

\ In spite of the progresses in non-perturbative techniques, the computation 

of hadronic matrix elements is still an open problem, particularly when the 
final state contains more than one meson. In this case, methods based on 



X 



Ph. Euclidean field theory, such as QCD sum rules or lattice QCD, have serious 

r-| ! difficulties in computing physical amplitudes [p], ||. Besides the standard 

parametrization of hadronic matrix elements in terms of B parameters, it is 
useful for phenomenological studies to introduce a different parametrization 
based on the contractions of quark fields inside the matrix element. In the 
following, we briefly discuss the Wick-contraction parametrization introduced 
within the framework of non-leptonic B decays in ref. ||. 

To be concrete, let us illustrate how this parametrization works in few 
examples taken from B physics. Consider the Cabibbo-allowed decay B + — > 
Z)°7r + . Only two operators of the AB = 1 effective weak Hamiltonian con- 
tribute to this amplitude, namely 

(D\ + \Qt c=1 \B + ) = (D\+\h,(l- l5 )duY(l- l5 )c\B + ), 
(D°<k + \Q£ c = 1 \B + ) = (D\+\h»(l-l 5 )cur(l- l5 )d\B + ). (1) 



*Talk given by M. Ciuchini at KAON '99, June 21-26, 1999, University of Chicago, 
Chicago, IL, USA. 
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DE(q 3 ,q 2 ,q 1 ;B,M 1 ,M 2 ) CE(q 3 ,q 2> q 1 ;B,M 1> M 2 ) CP(q 3 ,q 2> q 1 ;B,M 1> M 2 ) DP(q 3 ,q 2 ,q 1 ;B,M 1> M 2 ) 



Figure 1: Emission and penguin Wick- contraction topologies. The double 
dots denote the operator insertion. 



In this particularly simple example, the quark fields in the operators can be 
contracted only according to the emission topologies DE and CE, shown in 
fig. [I]. In the absence of a method for computing them, these contractions can 
be taken as complex parameters in phenomenological studies. The matrix 
elements can be rewritten as 

{D\ + \Q* C=1 \B + ) = CE LL (d,u,c;B + ,D°,n + ) 

+DE LL (c,u,d; B+, 7T+, D°) , 

(D°7r+\Qi c = 1 \B+) = CE LL (c,u,d;B + ,D°,7r+) (2) 

+DE LL (d, Uj c; B + , 7T+, D°) . 

The subscript LL refers to the Dirac structure of the inserted operators. In 
general there are 14 different topologies Of course, in order to be pre- 

dictive, one needs to introduce relations among different parameters given 
by dynamical assumptions based on flavour symmetries, chiral properties, 
heavy quark expansion, 1/N expansion, etc. This approach proves particu- 
larly useful for studying the AS = 1 B decays. For instance, let us consider 
the decay B + — > K + n® . Its amplitude receives contributions from all the 
operators of the AB = 1, AS" = 1 effective Hamiltonian. We consider only 
the matrix elements of operators which are both proportional to the largest 
Wilson coefficients C% and Gi and leading order in the Cabibbo angle. They 
read 

(K + it G \Q c 1 \B + ) = (K + n°\h,(l- l5 )sc^(l- l5 )c\B + ) 

= DP LL (c, Sj u;B + ,K + ,rt°), 
(K+n°\Q c 2 \B + ) = (ir + 7r |67 M (l-75)cc7^(l-75)s|5 + ) (3) 
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= CP LL (c } s,u;B + ,K + ,n°). 

The penguin contractions CP and DP are shown in fig. |l[ We stress the 
difference between penguin operators, which we have neglected here, and 
penguin contractions, which can contribute to the matrix elements of any 
operator. This kind of non-perturbative contributions, called "charming pen- 
guins" in refs. ||, could dominate AS = 1, AC = B decays because other 
contributions are either proportional to the small Wilson coefficients C3-C10 
or are doubly Cabibbo suppressed, as in the case of the factorizable emis- 
sion topologies of Qi2- A detailed analysis of non-leptonic B decays in this 
framework can be found in refs. |5], [7|. The presence of "charming pen- 
guin" contributions are likely to make the naive factorization approach fail 
in describing this class of decays. 

A different, but related parametrization of hadronic matrix elements has 
been recently proposed in ref. |J. In this approach, the parameters are the 
suitable combinations of Wick contractions and Wilson coefficients which 
are renormalization scale and scheme independent. In this way, the rela- 
tions among contractions enforced by the renormalization group equations 
are explicit. Besides, the phenomenological determination of the parameters 
do not depend on the choice of the Wilson coefficients. On the other hand, 
imposing relations among parameters based on dynamical assumptions may 
be more involved. 

Matrix-element parametrizations are less useful when applied to K de- 
cays, because there are few decay channels to fix the parameters and test 
the assumptions. 1 In addition, chiral relations allow the connection between 
matrix elements with one pion and those with two or more pions in the fi- 
nal states, the former being calculable with lattice QCD. However, a reliable 
lattice determination of (irir\Q 6 \K) , the dominant contribution to e'/e ||, is 
presently missing ||. 

Let us apply the Wick-contraction parametrization to K —>■ and verify 
whether there could be a connection between the longstanding problem of 
the AI =1/2 rule and a large value of the matrix element of Qq, as sug- 
gested by the recent measurement of e'/e. In terms of the Buras-Silvestrini 
parameters |J, the amplitudes K — > nn with definite isospin are 

ReA 2 ~ -(£ 1 + £ 2 ) ) 

1 Indeed, in the case of K — > 7T7t, there are only two complex amplitudes corresponding 
to the irir(I — 0, 2) final states. 
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ReA ~ (-^E 1 + ^E 2 -A 2 + P[ + P^j , (4) 
lmA (P x + P 3 ) , 

where E\^ 2 are the emission parameters, A 2 is built with annihilations and 

5 

L 

i=2 
10 



Pi = X!(?/2i-l(<52i-l)c£ + U2i(Q2i)DE 
i=2 

10 5 
+ J2(vi(Qi)cP + yi(Qi)Dp) +^2(y2i-l(Q2i-l)cA +U2i(Q2i)DA 



i=3 i=2 
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= X!(^2i-l(Q2i-l)/)A + y2i(Q2i)CA) (5) 
i=2 
10 

+ ^2(yi(Qi) cpa + yi{Qi)DPA) , 



i=3 



are the penguin-like parameters. The notation {Qz)ce refers to the connected 
emission with the insertion of the operator Qi, etc. 

Neglecting annihilations, as suggested by the large-iV counting or by 
CPS+chiral symmetries [|l(|, we are left with four parameters and three 
measured quantities. It is unlikely that ReA is dominated by emissions, 
since the large ratio ReA /ReA 2 would require large cancellations between 
Ei and E 2 , see eq. (f|). Therefore, in the most natural scenario, both ReAo 
and Irm4 are dominated by penguin parameters. Notice that Pi and P[ 
are different, so that no parametric relation between ReAo and Iim4 can be 
established. However, the following relation holds 

P[ = Zx(Qi)dp + z 2 (Q 2 ) CP + Pi(y -> z), (6) 

where jji and z t are the Wilson coefficients of the 3-flavour effective weak 
Hamiltionian. Given this relation, we can envisage a dynamical mechanism 
to connect the two parameters. Let us assume that P[ 3> Pi(y — > z) and 
that (Qi >2 )dp an d (Q$,6)dp share the same enhancement. Arguments may 
be provided to assume that [TIJ 

/ = (Qi)dp ~ l/N c (Q 2 ) CP ~ -<g 5 )up ~ -1/N c (Q 6 )cp ■ (7) 



2 We found that these two assumptions are compatible. 
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By using the experimental value of KeA and factorizing the emission con- 
tractions, we extract /, from which we derive 



B x = -9, B 2 = 7.5, B 5 = B 6 = 1.5. (8) 

It is interesting that the same mechanism enhances the B parameters entering 
ReA by a factor of ~ 10 and those of Imv4 by a factor of ~ 2, as required 
by the theoretical calculations to explain the experimental data. 

Alternatively, we could assume P[ ~ P\{y — > z) in eq. (§), namely that 
everything comes from the penguin operators Q§ and Qq. This is the old 



suggestion of S VZ Jl3| . Using perturbative coefficients, it is possible to show 
that this requires B s ~ 20 in order to fit (ReA )exp- Such a large value is 
excluded by the measurement of e'/e. 

To summarize, a connection between the enhancement of KeAo and a 
large value of e'/e cannot be established without some assumption on the 
long-distance dynamics. We have presented a simple example, which assume 
penguin-contraction dominance, that shows the correct pattern of enhance- 
ments. In this respect, models could give some insight, but quantitative 
predictions may prove hard to produce. Hopefully, non-perturbative renor- 
malization and new computing techniques will help overcoming the problems 
which prevent the lattice computation of Re^o and B Q [111] . 



M.C. and L.S. thank A. Buras for useful discussions and excellent steaks, 
beer and particularly desserts. G.M. looks forward to acknowledging the 
same in the future. 
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